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Abstract 

It is shown how the exchange coupling between two ferromagnetic planes 
embedded in an infinite non-magnetic metal, regarded as a function of the 
distance between the planes, may contain important components which oscil- 
late with periods not predicted by RKKY theory. The interesting case of a 

FCC(llO) structure with a Cu-like Fermi surface is discussed in detail. 
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Oscillatory exchange coupling between metallic magnetic layers across a non-magnetic 
spacer has been intensively studied over the last five years. One of the main issues in this 
area has been the determination of the oscillation periods of the coupling as a function of 
the spacer thickness. 

The physical mechanisms which have been proposed for explaining this phenomenon 
include the quantum well theory (QWT) of Edwards et al. [HHH an d an extension of the 
RKKY theory to the multilayer geometry due to Bruno and Chappert [fU|. In the QWT 
electrons propagating across the multilayer structure experience spin-dependent potential 
wells whose depths depend on the exchange interaction in the ferromagnetic layers. The os- 
cillatory behaviour of the interlayer coupling arises as a consequence of quantum interference 
effects inside the wells, and bears a formal analogy to de Haas- van Alphen oscillations [|TJ. 
The existence of these quantum wells has been confirmed experimentally by photoemission 
measurements ||. 

In the RKKY theory the oscillation periods are directly related to the spacer Fermi 
surface (FS) and are given by the wave vectors q z perpendicular to the layers that span the 
FS across those parts whose group velocities are mutually antiparallel. There is a general 
belief that all oscillation periods are given by the RKKY theory. Indeed, it has been shown 
that in certain simple models 0-01 the periods predicted by the QWT coincide with the 
RKKY ones and are given by the extremal dimensions of the spacer FS in the direction 
perpendicular to the layers. In another case JF|, where the lattice lacks reflection symmetry 
about a layer plane, the correspondence between the quantum well and the RKKY periods 
is more subtle but still obtains. Of course in the models mentioned above harmonics of the 
RKKY appear but no new fundamental periods. Furthermore, d'Albuquerque e Castro et 
al. [§] showed analytically for a very general model that RKKY theory holds in the limit 
of very small exchange splitting in the ferromagnetic material. However, van Schilfgaarde 
and Harrison found that real systems, such as Fe/Cr, are not in this limit ||, although they 
were convinced that the oscillation periods were derivable from the RKKY. Here we show, 
however, that under certain conditions, including those met in FCC (110) multilayers with 
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a Cu-like spacer FS, oscillation periods exist which are not predicted by RKKY. These are 
in addition to the usual RKKY periods. 

A FCC(llO) magnetic multilayer within the one-band nearest-neighbour tight-binding 
model is an interesting system because of the difficulty in determining analytically the 
energies of the resonances and size quantized states in such a structure [|10|| . Here, however, 



we approach the problem using the formalism of Ref . f| , which gives the coupling in terms of 
the one-electron propagators. This enables us to find the periods of oscillation analytically 
and to evaluate the exchange coupling numerically. 

We consider a multilayered system consisting of two parallel ferromagnetic atomic planes 
embedded in an infinite non-magnetic material. We label these two planes and n, so that 
the number of atomic planes in the spacer layers is equal to n — 1. It has been shown recently 
that as far as the interlayer coupling as a function of the spacer thickness is concerned, the 
thickness of the magnetic layers affects the phase and amplitude of the oscillations, but not 
the periods JTTJ . 



Within the single band model, the expression for the exchange coupling J, defined as 
the difference in the thermodynamical potential between the ferromagnetic and antiferro- 
magnetic configurations, is given by || 

J = --£ fduf(u)F( m ,u), (1) 

where 

F(q h u;)=Imln{l + W e lG{ (q h uj)GU<l\\,u)}, (2) 

and G^ (q||, is the off-diagonal matrix element between planes and n of the Green's 
function for an electron with spin a in the ferromagnetic configuration of the system, f(u) 
is the Fermi function, V ex is the exchange interaction in the ferromagnetic layers, and the 
summation over qn is restricted to the two-dimensional Brillouin zone (BZ). We assume for 
simplicity that the site energies in the ferromagnetic material are = e^V ex , where eo is 
the spacer on-site energy. In this situation, the off-diagonal propagators in Eq.(^) can be 
written for each qn and u as 



where t^i = iV^liV^goo) -1 ; and g n o is the matrix element of the bulk spacer Green's 
function. A similar expression can be obtained for Gq'^- In general, within the one-band 
model, g 0n is equal to g n0 , apart from a possible qy-dependent phase factor. 

It is clear that the behaviour of the coupling as a function of the spacer thickness is 
related to the dependence of g n o on n, which can be determined as follows. The matrix 
element of g between arbitrary planes I and m is given by 

gim = (;H / dq ± — — — (4) 

2n J-z- uj — E(q\\, q±) + 1 0+ 

where q± is the wave vector perpendicular to the layers, d is the interplane distance, and 
E(ci\\,q±) describes the bulk spacer band structure. We evaluate the above expression 
for gi m with / < m by integrating around the boundary of the semi-infinite rectangle 
—7i/d < Req± < ir/d ; Imq± > 0, in the complex g^-plane. This procedure is com- 
pletely general and can be applied, within the one-band model, to any lattice structure 
and layer orientation, with hoppings to arbitrary number of neighbours. For the FCC(llO) 
case, E(q\\,q±) = e(q\\) + 2t\(c[\\)cos(q±d) + 2t 2 (q||)cos(2gj_(i), where e(q[|) = —2t Q cos(2q x d), 
= —4tocos(q x d)cos(y/2q y d), and t 2 (q,\\) = —t . Here q x and q y are the components 
of qn, — to is the hopping between first nearest neighbour atoms, and the origin of energy 
is chosen such that eo = 0. ti and t 2 are the hoppings to first and second nearest planes, 
respectively. Using the contour integration described above, we obtain 

g n0 (qy ,u)=A 1 (qy , u)e<* ^ ^ nd + A 2 (q,, , u)e^ ^ )nd , (5) 

where cos{qjd) = + 8{u - e + 2t 2 )/2t 2 ]/4, Aj = [2i(cos(q 1 d)-cos(q 2 d))(l- 

cos 2 (qjd)) 1 / 2 ]~ 1 , for j = 1,2, and 7 = t\jt 2 . Here q± = ±q% and q± = ±q 2 are the roots of 
the equation E(q\\, q±) = uj. In the present case go n = g n o and the above expression for g n0 
can be extended to continuous values of n. Eq.([5]) shows that for values of qy and u, for 
which q 1 and q 2 are real, g n0 oscillates with the superposition of two periods, 2iid/\qi\ and 



2nd/\q2\, which are in general incommensurate. In those cases g n o exhibits a quasi-periodic 
dependence on n. As it is shown below, this fact may have a striking effect on the coupling. 

The function F in Eq.(|I]) varies with n through g n Q and therefore exhibits the same 
quasi-periodic behaviour. In order to deal with this quasi-periodic function, we make use of 
a procedure analogous to the one recently proposed (TIJ to investigate the dependence of the 
coupling on the magnetic layer thickness. It consists in replacing n in expression fl5|) for g n o 
by fictitious variables rii and ri2, which multiply q\ and q2, respectively. The real physical 
situation corresponds to n\ = ri2 = n. The extended function F(ni,ri2) is then periodic in 
each variable separately, and can be Fourier analysed in the usual way. Thus we find 

F= J2 C mum2 e l ^ +m ^ nd , (6) 

mi ,7712 

where C miiTn2 (qy, uj) are the Fourier coefficients, and m x and m 2 are integers. Since F is in 
fact a function of g^ , it follows that C mi ,m 2 — unless mi + m 2 is even. These coefficients 
contain all the information about the electron potential in the ferromagnetic layers. In 
particular, they depend on the magnitude of V ex and vanish for V ex = 0. In the RKKY limit, 
where F is replaced by the leading second-order term of its expansion in powers of V ex , the 
Fourier coefficients can be determined analytically. In this limit only six coefficients appear, 
namely C±2,o, Cq,±2, and C±± j ±i, whose values are given by 6*2,0 = V^ x Al/2i, Co,2 = Ka^l/^, 
C lfl = V^A^/i, and the property C_ mi _ m2 = C* mum2 . 

By inserting Eq.(|6]) into Eq.(jl]) we find that, according to the usual stationary phase 
method [p]-|^JTT[], for sufficiently large values of n, the nonzero contributions to the coupling 
come from u equal to the Fermi energy Ep and qii in the neighbourhood of those points at 
which the argument of the exponential is stationary. They are given by the equation 

m x V\\q^ h E F ) +m 2 V ll q 2 (c lll ,E F ) = 0, (7) 

where Vy is the two-dimensional gradient in qy space. The position of the extremal points 
can be determined exactly from the analytical expressions for q\ and q 2 . The weight of the 
contribution to the coupling from each solution of Eq.(|7|) depends on the magnitude of the 
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corresponding Fourier coefficient and partial derivatives of = ni\q\ + 7712^2 with respect to 
q x , q y , and u. 

Clearly the surfaces q± = ±qi(q\\, E F ) and q± = ±52 (qy, Ep) map out the FS. Fig. [I] 
shows a cross-section perpendicular to the atomic planes of the FS for E F /2to = 1.64 and 
q y = TT^/2/Ad. Note that q\ > and q2 < 0. Full lines correspond to the surface ±gi and 
dashed lines to ±q 2 . The nature of the FS depends on Ep, and three distinct energy regions 
are to be considered, namely — 12io < Ep < — Ato, — 4to < -E'f < 0, and < Ep < 4to- The 
two boundary values, — 4to and 0, correspond to those values of Ep at which the FS first 
touches the layer geometry BZ, and develops necks, respectively. As we show below, the 
numbers of periods we obtain in the three regions are different. 

It is interesting to examine first the predictions for the periods in the RKKY limit, where 
only six integers pairs m x m 2 are to be considered in Eq.(^). In the first energy region we 
find only one period A^q = 7rd/|gi(q|, E F )\, with qff = (0,0). In the second region an 
additional period = 27rd/|gx(q|, Ep) + q 2 (q^, E F )\ appears, where = (0, ±7r\/2/4(i). 
Finally, in the third region RKKY predicts four periods, namely, X^q, ^20 = ^d/lqi^q^, Ep)\, 
= nd/\q 2 {qf\,E F )\, and Af x = 27rd/| 9l (qjf, Ep) + & (q|, Ep)\ = \f A = 2nd/\ qi (q{ , E F ) + 
q 2 {q( l ,E F )\, with q| = (±7r/2d, ±7r^/4d) and qjf = (±7r/2d, 0) . The Fourier coefficients 
associated with the period A^ 1>m2 are C± mit ± m2 (c^, Ep). These periods are shown in Fig. || 
as functions of Ep. For the present model we find that ^(qy, Ep) and ^(qff, Ep), which are 
represented in Fig. [[], satisfy the relation gi(qjj, E F ) — ^(qy, Ep) = ir. Thus, the oscillation 
periods A20 and Aq 2 cannot be distinguished just by looking at discrete integer values of n. 
We recall that in FCC Cu E F lies in the third energy region with FS necks. A quantitative 
description of the oscillation periods for Cu can be obtained within the present framework 



by going beyond nearest neighbours and using the tight-binding parameters of Halse ||12 
Then A J 1 and \f 1 become distinct periods and, by taking into account the equivalence of 
A| o and Aq 2 for a discrete lattice, we find exactly the four RKKY periods of Bruno and 
Chappert [|]. 

However, it follows from Eq.(0) that there are contributions to the coupling with periods 
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other than those which arise in the RKKY limit discusssed above. In fact, it is easy to show 
that for Ep in the second and third energy regions we find an infinite number of solutions, 
corresponding to infinitely many values of m\ and m,2- Some of these solutions correspond 
merely to harmonics of the fundamental RKKY periods but some fundamentally new periods 
can arise. The simplest and most interesting case actually occurs in the top energy region, 
where Vygi and Vy^ vanish simultaneously at qjj. Thus, Eq.(^) is automatically satisfied for 
any values of mi and m.2- The corresponding periods are A^ m2 = 2 r Kdj\ra\q\ + m%q}\. Fig. 
H] exhibits the coupling J at temperature T = 2.0 x 10 W/ks as a function of the spacer 
thickness for Ep/2to = 1.64 and V ex = 0.15W, where W = 16£o is the spacer band-width. 
We chose this value of E F so that an important new period Ag l5 plotted as a function of 
Ep in Fig.|2|, is well separated from the RKKY periods. The full line in Fig.^| corresponds 
to the result obtained from Eq. ([!]), and the dashed line to the RKKY approximation scaled 
down by a factor of 8. In both cases n was treated as a continuous variable, but the physical 
discrete values are indicated. For this Ep the long period dominates both curves although 
with an amplitude differing by a large factor, but the interesting fine structure is different and 
reflects contributions beyond the fundamental RKKY periods. To make this point explicit, 
we present in the inset the absolute value of the ratio between some coefficients C mi , m2 and 
C\ t i, the largest coefficient of a fundamental RKKY period, as a function of V ex and for 
qjj. As expected, for very small exchange splittings, the magnitudes of higher order Fourier 
coefficients relative to that of the fundamental RKKY one are negligible. However, they 
rapidly increase with V ex , and the ratio |C mi)TO2 1 / 1 Cx,i | becomes significant. From the inset 
in Fig.0, we see that for V ex = 0.15W, there are additional contributions to the coupling 
coming not just from harmonics of the RKKY frequencies, which correspond to m\ and 
m 2 even, but a very important one coming from a new period A3 1 — 4.5d. In fact, the 
contribution from this new period can be calculated separately using the stationary phase 
method p|-^,pT|1. The result is shown in Fig.|], together with those corresponding to two of 
the fundamental RKKY periods, namely, \% and A^q. It is clear that this new period is as 
important as the RKKY ones, except for the dominant long period whose large amplitude is 
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due to FS geometry which causes some second derivatives of to vanish, making a stationary 
phase evaluation of this contribution impossible. The amplitude of the new contribution falls 
off as 1/n 2 just like the normal RKKY components and their harmonics [0-0. 

The appearance of the non-RKKY periods is clearly related to the spacer off-diagonal 
propagators oscillating as a function of the spacer thickness with more than one period, for 
fixed energy and qn . As we have shown, such a behaviour can be found even in the one-band 
model, for which the FS is simple and has a single sheet. It can also be found in those cases 
in which the spacer FS has more than one sheet. Thus we may expect the occurrence of 
non-RKKY periods in the coupling through a non-magnetic transition metal. Therefore, 
in those cases, the interpretation of the results in terms of just the RKKY theory may be 
misleading. 

In conclusion, we have shown that the relation between the oscillation periods of the 
coupling and the spacer FS is more complex and subtler than has been assumed so far, 
making possible the appearance of non-RKKY periods. This is the central result in this 
communication, which settles the long standing question about whether or not the quantum 
well and the RKKY theories always give the same periods of oscillation. 

We are grateful to EPSRC and Royal Society of UK, and CNPq of Brazil for financial 
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FIGURES 

FIG. 1. Cross-section of the spacer Fermi surface perpendicular to the layers for Ep/2to = 1.64 
(see text). Full lines correspond to the surface ±gi and dashed lines to ±q2- The two vectors 
qi(q^,Ep) and </2(qfj, -ETf) are indicated. 

FIG. 2. RKKY periods as a function of Ep. Curves 1,2,3,4, and 5 correspond to Af^, l5 
Ago, Ag 2 , and Af ;1 = \f 1} respectively. The dot-dashed curve is a new period \\ x . 

FIG. 3. Exchange coupling as a function of the spacer thickness for Ep/2to = 1.64, V ex = 0.15W 
and fcfiT = 2.0 x 1CT 3 IF (full line). The dashed line corresponds to the RKKY result scaled down 
by a factor of 8. The inset shows the ratio |C 3jl | /| Cx,x | (full line), | C*4,2 1 /I Ci,i I (dashed line), and 
|C4,o|/|Ci,i| (dot-dashed lined) as a function of V ex , for q|j. 

FIG. 4. Contributions to the coupling coming from the non-RKKY period A3 1 (full line), and 
two RKKY periods Af (dashed line) and A^o (dot-dashed line) (see text). 
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